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ABSTRACT 

This evaluates the United States Air Force (USAF) Nonnuclear Armament 
Program (NAP) models and specifically the Heavy Attack (HA) model. Particular 
attention is paid to the optimization techniques incorporated in Heavy 
Attack, to the validity of the inputs being optimized, and to the implications 
of underlying model assumptions. An examination is made of the validity of 
using target values as model inputs for not only the beginning of a conflict, 
but also for times extending into the conduct of a conflict. New technoloqy 
has been applied to the problem and the success achieved is reviewed. 
Reformulations aimed at improving model capabilities and/or solution speeds 
are described. 
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I. THE PROBLEM 



Military Planning for the United States primarily involves preoaration 
for an aggressive response to an attack made against the United States or 
allied forces. This basically defensive posture all but denies the 
United States the advantages inherent in surprise, and requires, as a 
counter to the initiative granted an attacker, that the United States 
maintain operational effectiveness and vigilance second to none. 

To achieve maximum military effectiveness given fiscal constraints 
requires careful consideration in planning and budgeting. Men without 
weapons will neither intimidate nor stop an enemy, neither will men with 
weapons but without munitions or a means of maneuver (transport). As an 
example, the United States Air Force must weigh each dollar spent to see 
if it should be best utilized for personnel (acquisition, retention, and 
training), aircraft (acquisition and maintenance), fuel, or munitions in 
order to best fulfill its many missions. 

The research described herein addresses the Air Force nonnuclear 
munitions procurement model. This model requires as input: 

1. A target list to include target military values, defensive capabil- 
ities, and factors describing the potential for confirming target 
kill; 

2. Numbers of friendly aircraft and munitions available; 

3. Probabilities of encountering differing weather conditions; and 

4. The effectiveness of specific munitions against various targets 
when delivered in a weather condition by a specified aircraft. 

Given these inputs, the model seeks to maximize the damage done to 
the enemy by planning the use of the most effective munitions. 
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The model is played separately for each major Air Force theatre. 
Individual theatre results comprise the single largest input for non- 
nuclear weapons procurement decision-making for the Air Force. The 
amount of money spent with the aid of this model is currently in excess 
of one billion dollars a year. Not only is the model central to Air 
Force budgeting and planning, but it is also relied upon by the major theatre 
commanders for insight into scenarios involving their current missions. 

With so much at stake in its use, an examination is in order to 
better understand the assumptions, formulation, processing, mathematical 
solution, and solution report generation for this model [Ref. 1]. 
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II. DESCRIPTION OF THE HEAVY ATTACK (HA) MODEL 



Heavy Attack (HA) is a program which uses an imbedded nonlinear 
optimizer to identify a set of sortie allocations which maximizes the 
military worth of targets killed. It has been used since 1974 as part of 
a set of computer programs that have been known as Saber Mix or, more 
recently, as the Nonnuclear Armament Plan (NAP) models. These models 
together attempt to provide an optimal munitions mix for a given 
specification of available sorties, targets, and other factors. 

The HA model can be viewed as consisting of an internal (optimization) 
model and an external model consisting of input and the output sequences. 
This view will be adopted in describing the model and its variations 
primarily because the optimization is still modeled, if not processed, 
exactly as it first was eight years ago, while input and output models 
have changed (grown) continuously. Careful attention will be paid to the 
optimization model while the input and output models will be examined 
only to achieve an understanding of the consequences of model input 
aggregation and output unraveling for solution interpretation. 

A. DESCRIPTION OF THE HA OPTIMIZATION FORMULATION (INTERNAL MODEL) 

HA solves a sequence of internal models. Each of these models has 
the same mathematical structure (formulation), and each is a nonlinear 
optimization problem by virtue of its objective function. The internal 
optimization model is formulated as follows: 
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maximize: 
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- Dj) 




J 




subject to: 


0?1 


Sij = s,, t = 1, 






KjlTj, j = 1, 



(E (1 - V Si j + E (-V s, j)<0, 

\J'0„ <«■’ <«■>/ 



hi > «• 



m = 1, 2, . . . M; 



The objective function, which is convex with respect to S.^ 

[Ref. 2: pp. 8, 9], quantifies the value of all targets killed in a 

time period as follows: 

E Vj(Kj<Sij> - Dj). 



where, 

i 

j 






V . 
J 



= aircraft type index (i =1, ..., I); 

= target type index (j = 1, ..., J); 

= number of sorties, the independent vaviables; 

- value (military worth) of target type j; 

= cumulative number of targets killed in prior time periods; and 



= number of kills of target type j, a nonlinear function of S- 

J J 






K- is defined to be a function of the independent variables S--: 

J ' J 

= ^(1 - exp ( - + Z ?iiS, .,))), 

J 



Tj ' 0 ft -'no' 



where. 



T. = number of type j targets; 

J 
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C. = target kill confirmability parameter (0 _< C- £ 1) controlling the 
^ extent to which the law of diminishing prodactivity (as described 
later) applies; 

T. C. 

0. = ^ log (1 - yi 0.) ; a term (to be explained) added for mathematical 

^ j ^ convenience; and 

P-- = expected number of type j targets killed per type i aircraft 
^ sortie when no other targets of type j have been previously 
killed and when conditions of kill confirmability are perfect. 

HA problem constraints are of three forms: 

1 . Sorties available : 

J 

S.. = S^., i = 1, 2, I; 



where S- = the number of sorties for aircraft type i. 

2 . Target : 

1 1 Tj, J = 1, 2, J; 

where = the lower bound on targets of tyoe j which must be 
killed.^ 

3. Flight composition with the general form: 

^m ( " ®m^ ^i j ^“®m^ ^i 

VjcJ^ mJ/ 

where, 

m = flight composition constraint index (m = 1, 2, M); 

i|^ = aircraft type; 

'5^ = +1, maximum (-1, minimum); 

= set of targets for which a maximum (minimum) flight 
composition is required; and 

= maximum (minimum) proportion of sorties flown by aircraft 
type i against targets included in set J 
[Ref. 7: pp. 5, 8, 101. ^ 

Kj is computed assuming "a law of diminishing marginal productivity . . . 

i.e., the number of targets per sortie decreases for each successive sortie" 
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[Ref. 2: p. 1]. This underlying assumption is based on a belief that 
targets will be harder to find, harder to surprise, and harder to kill as 
the battle continues. 

By restating target constraint inequalities to solve for z 
the problem can be formulated (see Appendix A) with linear constraints 
[Ref. 2: pp. 6-7, 10-12] of the form: 



ij< E 

J i=l 



<«J 



where, 

-T. 

-3 " ~r 

iTj (1 - ‘j tf >> - ”j- 

J J 

Further examination of K- provides significant insight into the 

xJ 

nature of the function which is being maximized. If 

T. C, I 



K 



. .^(1 . exp . E PijS,j))). 



then noting, 

exp = 1 - 

J J 

and (for convenience of exposition) letting. 



the term - Oj (from the objective function) can be restated (see 
Appendix 8) as follows: 

T. C. 

Kj - Oj = ( ji - Dj)(l - exp (XJ (- J )))• 

J J 
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The function K- - 0^ is a composition of two embedded functions. 

J J 

The constant parameter C- can be used to select either component. When 
Cj is equal to one, 

K. - Dj. = (Tj - 0j)(l - exp ^ ). 

When C. is zero. 




Regardless of C., if X./T. is small (a target-rich environment exists), 

J J J 

then K . - 0 . = X . . 

J J J 

C. can be interpreted to be a coefficient of confirmability. That 

J 

is, if conditions allow a pilot to confirm the effect of his first 
ordnance drop before delivering any others then C. should be zero. This 

J 

case can be called shoot-look-shoot. 

If conditions are such that a pilot must deliver his ordnance without 
reference to the success or failure of any weapons that might have been 
delivered earlier, then should be one. This case could be called 
dump-all-ordnance. 

While the composition of K- - 0- can be understood for the boundary 

J tJ 

conditions where C- is equal to one or zero, it is not at all certain how 

J 

K. - D- should be interpreted in cases where C- lies between the two 

J J J 

boundaries, except that it will behave as some mixture of the two embedded 
functions. 

Summarizing, HA solves a problem with a nonlinear, but convex, 
objective function subject to three sets of linear constraints. Function- 
ally it selects sorties, not weapons j to inflict the greatest damage upon 
an enemy. Constraints demand each aircraft be utilized when available, 
that a specified range of each target type must be killed, and that 
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specified sets of targets, associated with an aircraft, must be attacked by a 
minimum (or maximum) percentage of the sorties available for that aircraft. 

B. DESCRIPTION OF THE JUNE 1982 HA EXTERNAL MODEL 

Because the optimization model is formulated in terms of sorties 
( aircraft/target combinations) while the purpose of HA is to provide 
preferred weapon information, the HA external model largely concerns 
itself with building sortie information (P^j's) from inputs and extracting 
from sortie solutions (S--'s) information regarding preferred weapons. 

The HA model, external to the optimization model, has been modified 
repeatedly in the years since 1974 for a variety of reasons. Determining 
the history of these changes might prove interesting, but what is more 
important is to understand HA in its present form. (Inputs, processes, 
and outputs for the original HA model, as described in the seminal paper 
for the model [Ref. 2: p. 2], are described in Appendix C.) 

Preceding HA in the execution of the Nonnuclear Armament Plan (NAP) 
model programs are Weaponeer, Survivor, and Selector [Refs. 3: p. 5 and 

Ref. 4: p. 2]. 

"The first model. Weaponeer, computes the expected target kills 
per pass for various aircraft/munition/target combinations. . . . 

The second model. Survivor, computes the attrition of delivery 
aircraft on an iterative basis and determines the expected kills 
per sortie .... The third model. Selector, selects the preferred 
long-list weapon for each aircraft/target/del ivery band [weather 
type] combination. ..." [Ref. 4: p. 2] 

The models preceding HA in the NAP sequence are of current interest 
only insofar as they provide inputs to the HA optimization. What should 
be noted, however, is the concern for the passage of time implicit in the 
short descriptions of Weaponeer and Survivor given above (and in the 
prior discussion of diminishing marginal productivity), as well as the use 
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of weather types in Selector. While Clasen, Graves, and Lu [Ref. 2] 
never discuss use of the model over multiple time periods or delivery 
(weather) types, these two elements are currently incorporated in the HA 
model. (Multiple time period use is implied in [Ref. 2] by inclusion of 
the term D- in the objective function and by the adjective "interval" 

J 

describing model inputs and outputs, however, no "iteration-over-time" 
scheme is described. The effects of weather on sortie effectiveness are 
never mentioned in [Ref. 2].) 

Weather conditions have a profound impact on aircraft/weapon effective- 
ness. The HA external model expresses weather condition as six discrete 
"bands," associated with extremely poor to essentially unlimited flying 
visibility. The June 1982 HA input model builds sorties and associated 
expected kills using a weighted average based on expectation of weather 
types in the theatre being modeled. The P. • associated with each S- • is 

* J * J 

formed as follows: 



May 

^ Pw T Pijkw 



Mumv 
’J W=1 



where, 

k = set of all weapon types; 

p^ = probability of a particular weather type; 

P^jkw " expected number of kills by aircraft i loaded with weapon k 
against target j in weather condition w; and 

numw = number of weather condition types modeled. 

The effect of this averaging prior to the optimization model is two- 
fold. Sorties are not limited to a single weapon type; that is, sorties 
are evaluated as aircraft, not as aircraft/weapon combinations. More 
importantly, an average quantity representing the aircraft's effectiveness 
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is being sent to the optimizer. The result is a selection of extremal 
averages rather than the averaging of selected extremal values. 

This can lead to paradoxical recommendations. Consider a paradigm 
with two equally probable weather conditions, two aircraft types, and a 
single target type. Given one aircraft is moderately effective against 
the target in both weather types while the second aircraft is near 
certain to kill the target in one weather type and as certain to miss it 
in the other, the first aircraft is likely to be selected by the model 
(to be flown in both weather types) despite the dominance of the other 
aircraft in the second weather type. 

In HA the effects of weighted averaging may be subtle. A situation 
is likely in which for four of six weather types, the effectiveness of an 
aircraft/best weapon combination against a target type is quite high 
while in the other two weather types, the effectiveness of this aircraft 
is negligible. If another aircraft has mediocre success in all weather 
types against the same target, it might well be selected by the optimizer 
because the optimizer is not allowed to see the "extreme" effectiveness 
of the first aircraft against the target type in particular weather 
conditions. A weighted average model will not necessarily provide valid 
answers for any HA scenario. 

A second major function of the external model is to accommodate 
processing over multiple time periods. The number of time periods is 
limited to seven in the June 1982 model, and each time period is typically 
defined to last fourteen to thirty days. Model parameters do change over 
time. The number of dead targets (0,.) is accumulated following optimization 
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for each time period. Other parameters (in particular, target values) 
are reinitialized at specified time intervals. 

The use of multiple time periods in HA is myopic. That is, optimiza- 
tion proceeds forward in time, in one pass, with no backtracking. (Each 
optimization seeks to inflict maximum damage upon enemy targets in that 
time period without regard for the effect "decisions" made in the current 
time period may have on outcomes for future time periods.) 

It is apparent that a target type with greater than average defenses 
will generally have a higher than average value. Having a higher target 
value will increase the likelihood the target will be attacked--at 
least in a short duration "conflict" — since we are simply trying to 
inflict the greatest damage upon the enemy. However, by avoiding this 
target type in a longer duration conflict and killing other less valuable 
and less heavily defended targets, the attacking aircraft might be used 
to inflict more damage on the enemy, over time, by surviving longer and 
killing a large number of less valuable targets. 

Lack of logistical constraints poses the most significant difficulty 
of a practical nature for the June 1982 model. HA uses a preferred 
weapon without regard to its actual availability. The internal model 
never "sees" weapons — it only "sees" sorties. During a time interval, 
the model will continue to "fly" a sortie type as long as the objective 
function value is improved (subject to the model constraints). Even when 
it is evident the supply of a preferred weapon has been exhausted in 
previous time periods, there is no way to prevent the June 1982 model 
from continuing to use that preferred weapon in subsequent solutions. 
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C. DESCRIPTION OF THE JUNE 1982 HA PROGRAM CODE 

Appendix D is a hierarchy chart of the June 1982 HA programs. Sub- 
program NONLIN and all those programs strictly subordinate to it on that 
chart comprise the nonlinear optimizer, while the other programs manage 
input and/or output. 

Calling disciplines, variable naming conventions, and documentation of 
the HA code are neither standardized nor consistent. Variable names in 
the system frequently do not agree with current system documentation or 
with [Ref. 2], nor are they consistent among subroutines. (Variable 
names within the optimization code are, for the most part, consistent 
with [Ref. 2].) Appendix E details the name changes of some of the 
variables pertinent to the interface between the main program and the 
optimizer. Appendix F [Ref. 4: pp. 4-5] provides a data flow for the 

entire NAP process as it was executed in June 1982. 
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III. DISCUSSION OF TARGET VALUES IN HEAVY ATTACK 



A. CRITICISM OF THE USE OF TARGET VALUES 

In a recent analysis of HA by G. Jenkins [Ref. 5] (an Air Force 

civilian employee and user of HA) the manner in which inputs are prepared 

for HA, as well as the optimization technique used in the model, are 

discussed. He concludes the optimization is straightforward and infers 

all inputs except one are derived objectively and correctly. The one 

input he expresses concern over is target values. 

". . . the entire process is based on optimization of firepower 
scores. This model calls it military worth which is probably more 
correct, since command centers and runways really don't fit in the 
context of firepower scores as do tanks. Nevertheless the purpose 
of this entire process, which is pursued with meticulous objectivity 
throughout each set . . . boils down to optimizing the relative 
subjective worth of target values. Granted, this methodology is 
probably more credible than the proverbial smoke-f il led-room 
approach; however, it indicates that there is still room for 
improvement" [Ref. 5: p. 12] 

He continues, 

". . . there is room for concern over whether the model's solution 
is credible at all. Maybe the model merely serves as a guide to 
quantify and substantiate some decision-maker's intuitive feeling, 
so that he may proceed with what he always wanted or 'knew' to be 
true. Or, maybe there is a sincere interest in gaining insights 
into the combat process. If so, there is plenty of room for 
improvement in this model. . . . "[Ref. 5: p. 13] 

Jenkins is not alone in expressing concern over the use of target 
values in computing an optimal munitions mix for air forces. 

NATO first implemented the entire Sabre Mix (now called MAP) method- 
ology with few alterations, but later removed target values from the 
formulation. The new formulation simply seeks to maximize the total 
number ot targets killed subject to a (new) constraint for each target 
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type mandating predetermined target type kill proportions relative to 
other target types (a proportionality constraint). The new objective 
function, new constraints, and other problem constraints (much as 
originally found in HA) are optimized with a linear programming code 
[Ref. 6: pp. 4-11]. The reasons for embarking on the reformulation were 

given as follows: 

". . . While it is recognized that the concept of military worth 
plays a central role in tactical air mission planning models with 
short planning periods, this approach was abandoned for the following 
reasons : 

1. Considerable difficulty arose from trying to assign credible 
military worth functions to the different target types. 

2. It was considered advantageous to let proportions in which 

targets of different target types are killed constitute input 
to, and not output from, a model designed to contribute to 
the solution of a logistical problem." [Ref. 6: p. 3] 

The concerns expressed regarding the use of target values in HA can 
be summarized as follows. Target values are subjectively derived and, 
perhaps, invite manipulation of the model by those executing it. They 
are hard to assign in a credible fashion. They remain the only subjective 
input in an otherwise objective process, and their use might confuse the 
distinction between tactical and logistical decision making. 

Discussions with Major F. Cooper, the officer currently charged 
with executing the NAP models [Ref. 7], reveals another perspective on 
the use of target values in HA. The HA model is run separately for each 
major theatre (Europe, Pacific, S.E. Asia, etc.) by a theatre project 
team consisting of modeling personnel (from Major Cooper's office, 

AF/XOX) and military contingency planners and intelligence experts 
currently assigned to the command responsible for the theatre being 
examined. Projected targets are identified and target values are 
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assigned for the first time in a given scenario. Targets are evaluated 
using a scale of zero to twenty (relative to a tank platoon which is 
given a base score of one). The model is then executed for one time 
period. Based on the results for that first time period and an examina- 
tion of the projected replenishment capabilities of the allied and 
opposing forces, target values are assigned for the second time period. 

This process is repeated until the model has been run for all time 
periods. This cycle — process, evaluate, modify, and reprocess — has 
historically developed credibility for the model among its theatre users. 
Manual intervention also allows target values to be manipulated so as to 
discourage inadmissable solutions produced by the June 1982 HA model 
(because of missing constraints). Such manual, judgmental manipulation 
is done with the full knowledge and concurrence of the theatre project 
team, and its impact on model realism is carefully evaluated before 
continuing with the next time period. This careful evaluation of results 
for each time period often requires twenty or more model runs for a given 
scenario. However, the model is thus not allowed to use weapons no longer 
procurable, or to use weapons in quantities greater than can be procured, 
or to fly an aircraft to targets outside the aircraft's range. When a 
theatre study is completed, the results reflect the project management 
team's consensus and its total combined military judgment. 

8. A CASE IN BEHALF OF THE USE OF TARGET VALUES IN HA 

The use of target values in HA is viewed in vastly different ways 
by knowledgeable people familiar with the model. In reviewing the 
methodologies used by USAF nonnuclear munitions mix models, Jenkins found 
the use of target values in HA to be a weak link in an otherwise strong 
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chain. Loritzen finds their use difficult and inappropriate, while 
Cooper finds them absolutely necessary if the model is to reflect both 
the commander's priorities (as understood by his contingency planners and 
intelligence experts) and real-world constraints not accommodated by the 
model's formulation [Refs. 5, 6, 7]. 

Jenkins infers HA should replace the use of target values with a 

two-sided game reflecting optimal strategies, 

". . . we can see that this model is not totally unlike the other 
models .... However, there is no direct gaming structure. This 
model is purely a one-sided affair and the only service provided by 
the opponent is the supply of targets and an unaffected attrition 
rate. The staging of the scenario and reconstitution features of 
the model could involve certain optimal strategy games, but there 
is no evidence of this in the process . . ." [Ref. 5: p. 12] 

While the above is intriguing, developing the two-sided game preferred 
by Jenkins would require a scale, a set of values with which both sides 
will be originally endowed and by which each side will either gain or 
lose, depending on their strategies and initial endowments. It is 
apparent that this scale of values will have to relate to target values. 

If it is granted that values for potential U.S. targets are not derived 
easily, then how much harder must it be to evaluate targets on a scale 
applying to both the U.S. and its opponent(s)? 

NATO removed target values from the formulation for its strategic 
aircraft munitions mix model, but then apparently had a difficult time 
validating the proportionality constraint used in the new formulation. 
Studying the linear programming dual of their formulation reveals the 
dual variables are "marginal implied military worths of the targets" 

[Ref. 6: p. 15J. In fact, Loritzen concludes that in order to get valid 

results from the new formulation "the adjustment of proportionality 
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factors will necessarily have to be done manually by military experts 
until a situation is reached where the implied military worths do not 
differ significantly from their estimates ..." [Ref. 6: p. 15]. 

If reformulating HA as a game or reformulating it using target 
proportionality constraints provides no relief from the use of target 

V 

values, then perhaps running the engagement as a simulation model, from 
the first day of the engagement, can provide a more objective model input. 
It is this author's opinion that while the use of simulation model outputs 
for optimization inputs may be more objective, it is doubtful anyone 
would vouch for the output produced. Further, the cost of running such a 
simulation in all its required iterations might be prohibitive. 

Finally, an analytical model might be used to provide input in lieu 

of "subjective" target values; however, some significant, though more 

subtle, aspects of battle would prove particularly difficult to model. 

In evaluating alternatives for incorporating target activated munitions 

in the NAP models, Cudney and Bloomquist state, 

"The development of value curves must be based upon the judgment 
of military commanders and analysts who are experiened in esti- 
mating the effects that casualties and delay might, have on the 
success or failure of specific military missions." [Ref. 8: p. 6 

Combatant morale is also directly related to choices regarding target 
destruction priorities. Military commanders and analysts are the sole 
credible source for evaluating the effects of such factors. 

Without consideration of such factors as delay and combatant morale, 
this author would find it hard to accept any substitute for target 
values produced by an analytic or simulation model. 

^Emphasis added. 
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If no apparent, preferable alternative to the use of target values 
exists, why is their use viewed so negatively? It is likely that the 
practice of manipulating target values to compensate for model short- 
comings contributes to their poor reputation. If target values represent 
just that, the value of a target, then military analysts would find their 
use easier to accept as valid. It is difficult to stand behind model 
output when inputs must be manipulated so violently. 

Still, the issue of subjective model input remains. It is this 
author's belief the use of such input is allowable, even necessitated, 
because the use of military judgment is still the preferred alternative 
when trying to account for all the imponderables existing on a battlefield. 
Unfortunately, while admitting military science is a "soft science," 
analysts still seek to "exorcise" their models of subjective military 
inputs (judgments) rather than incorporating and exploiting their use. 
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IV. HEAVY ATTACK OPTIMIZATION ALGORITHMS AND THEIR PERFORMANCE 



A. THE ORIGINAL (JUNE 1982) NONLINEAR PROGRAMMING ALGORITHM 

The optimization code present in the June 1982 HA program (consisting 

of those programs strictly subordinate to subprogram OPTMUM in Appendix 0) 

appear to be a faithful translation from the original FORTRAN (to PLl 

and back to FORTRAN) of the optimizing code described in [Ref. 2]. That 

code solves a general nonlinear program by a sequence of local linear 

programs. The algorithm employed, 

". . . is a 'local,' 'gradient,' 'stepwise' correction descent 
algorithm .... By a 'stepwise' procedure we mean that given a 
point y^ in the domain of the functions, a 'correction' vector ay 
is determined and a new point y = y^ + k ay is used for the suc- 
cessor 'step.' It is a 'local' method because the correction 
direction ay and its length (determined by the skalar k) are 
obtained from the behavior of the system in a 'sufficiently' small 
neighborhood of the current point y . It is a 'gradient' technique 
inasmuch as the gradients of the function gi(y) are principally 
used to obtain the correction direction. . . ." [Ref. 2: pp. 13-14] 

Use of this algorithm to perform an HA optimization process, while 
precise and correct, involves significant expenditure of computing 
resources (typically five to six CPU minutes on an IBM 3022). Anticipat- 
ing a state-of-the-art optimizer would perform the optimization at less 
expense, the Air Force requested [Ref. 1] the June 1982 HA optimizer be 
replaced with the X-System. 



B. THE X-SYSTEM 

Like the original optimizer, the X-System solves a nonlinear problem 
with a sequence of local linear programs. (The HA internal model consists 

p 

The objective function. 
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of a nonlinear objective function subject to entirely linear constraints 
so that for each linear program, all that is required is that a linear 
approximation be made of the objective function and then a "standard" 
linear program is run.) 

The X-System has been operational, but under continuing development 
since 1974 [Ref. 9]. It is a general-purpose, state-of-the-art optimiza- 
tion system which is used as both a vehicle for research and as the basis 
for a number of commercially installed, customized applications optimizers. 
It consists of open FORTRAN subroutines and is implemented in FORTRAN IV. 
The subset of FORTRAN with which it is coded is accepted by a majority of 
FORTRAN compilers. 

The X-System is designed to solve large-scale optimization problems, 
and is especially effective on mixed integer problems. Decomposition 
issues have been an area of major interest to the designers; however, the 
core linear programming module has received the most design effort. It 
exhibits many unique features including: 

1. Hyper-sparse data representation [Ref. 10]; 

2. Complete, constructive degeneracy resolution [Ref. 11]; 

3. Basis factorization [Ref. 12]; and 

4. Elastic range constraints [Ref. 9]. 

In order to best support the wide variety of applications using the 
X-System, the system is designed to support all other optimization features 
simultaneously with the nonlinear feature (e.g., sortie constraints form 
an intrinsic Generalized Upper Bound, or GUB, set which has been exploited 
in the nonlinear solution). 
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C. NUMERICAL EXPERIENCE WITH THE X-SYSTEM IN HA 



The X-System was imbedded in HA as a subprogram, without changing 
the internal or external models, and delivered to the Air Force in 
July 1982. Headquarters, United States Air Force (XOX/FM) has been 
testing this version of the HA code {called "Fast" Attack by its users) 
since then. The X-System returns solutions whose objective function 
values agree with those provided by the June 1982 code to the second or 
third significant digit in less than one fiftieth of the time while using 
a fraction of the compute region. However, the numbers of particular 
sorties chosen differ significantly in some instances. These differences 
result from the level of precision specified for the optimal objective 
function value which permits early termination of the algorithm with an 
acceptable solution. 

Since both the optimizers described are supposed to deliver correct 
solutions, one can hardly help but be surprised when their optimal solu- 
tions for the same problem differ, no matter how small the difference. 

One fundamental difference in the two algorithms is that while the 
original optimizer included a coded gradient function for the objective 
function, the X-System uses an automatic numerical difference approxima- 
tion to estimate gradients. Use of this approximation in the X-System to 
enhance robustness, that is, to eliminate the errors and frustration 
associated with coding derivative functions, appears to be responsible for 
the small differences in objective function value. 

A solution precision factor is used by the X-System. This algorithm 
parameter directs the X-System to stop optimization at the first point 
where the tolerance of the solution is estimated to be comparable to user 
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requirements, or confidence in the precision of input data. (This avoids 
wasting computing resources extracting the "last few pennies" in optimality 
from a problem for which input coefficients have been rounded to the 
next higher dollar.) By tightening this factor somewhat, perceived 
instability in total numbers of individual sorties can be eliminated. 
However, if this tolerance factor is out of proportion with user confidence 
in input data, any "stability" achieved is illusory. 

The X-System has enabled HA to be used on a timesharing system. 
Internal models with 81 constraints (13 GLIB) and 793 variables typically 
yield solutions in 10-15 CPU seconds on an IBM 3033 processor, using 
approximately 250K bytes. The much faster response of the enhanced 
system has not only enabled speedier evaluations for the theatres, but 
has encouraged a critical review of the HA model. It is now technically 
possible to add model enhancements because neither space nor time constrain 
the problem as they did with the prior optimizer. 
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V. A NEW HEAVY ATTACK EXTERNAL MODEL (WEATHER WARS) 



In an attempt to ameliorate some of the major problems in the HA 
model, a new external process was designed in July and partially imple- 
mented in August 1982. It is still undergoing testing and refinement. 

The new model, known as "Weather Wars" (WW), does away with the use 
of weighted average sorties and attempts to logistical ly constrain the 
problem. 

WW avoids the use of weighted averages by building sorties for a 
particular weather type. Optimization is performed and inputs are 
updated and reinitialized as in the former HA; however, at the end of 
processing for the last time period, another "weather war" is "fought" 
in which sorties are built for a different weather type. 

Running separate "wars" for each anticipated weather type (typically 
six) is costly. However, the expense is offset by the speed with which 
the X-System provides internal model solutions. When weighted averaging 
of weapons expended is performed following completion of all "weather 
war" processing, the resulting weighted average solution is provided and, 
at no additional cost, so is the maximax solution, the number of weapons 
required to meet all constraints and inflict great damage upon the enemy 
even if the weather becomes the enemy's consistent ally. 

WW logistical ly constrains the problem by noting maximum procurable 
quantities for each weapon type and by reading in a user-defined number 
of "best" weapon types for an aircraft operating against a oarticular 
target type in a given weather condition. Prior to each optimization. 
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sorties are built and probabilities of effectiveness (P^^'s) are 
assigned using the "best" weapon type still available. If all the "best" 
weapon types for a particular aircraft/target combination are no longer 
available, the model builds sorties using a fictitious weapon type with a 
zero coefficient of effectiveness. 

The internal model still has no explicit logistical constraint, and 
given time periods of fourteen to thirty days, will include sortie totals 
in some time period solutions which consume weapons in quantities greater 
than will be available. However, this "over-use" of a weapon type can be 
made arbitrarily small by reducing time period duration (and subsequently 
increasing total numbers of optimizations). Another possible method for 
strictly limiting the numbers of a weapon type used is to run the entire 
WW model iteratively, reducing the maximum procurable quantity of a 
weapon type to the number used in the time period prior to the one in 
which the "over-use" occurred. 

WW permits analysis of realistic scenarios in which sorties are to 
be selected after weather conditions are known. The June 1982 HA infers 
an assumption that sorties must be planned with only synoptic weather 
forecasts, or that various weather conditions will exist throughout the 
theatre in specified proportions. WW presents opportunities for decision 
analysis incorporating weather-dependent target values and damage 
confirmabil ity. 
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VI. POTENTIAL FUTURE DIRECTIONS FOR HEAVY ATTACK 



A. HA INTERNAL MODEL AS A GENERALIZED NETWORK 

It is often advantageous to reformulate an optimization problem into 
an equivalent model which may be easier to solve. Recent computational 
advances in the efficient solution of generalized networks [Refs. 13, 14] 
have allowed these relatively specialized linear programs to be solved in 
a fraction of the time required to solve them with "standard" linear 
programming techniques. 

The internal HA model can be viewed as a generalized network if 
flight composition constraints are ignored and the objective function is 
simplified (by letting C. = 0 for all j, or assuming a target rich 

J 

environment) as follows: 

J I 

maximize:. ^ 

j=l i=l J 'J 'J 

The model can then be viewed as displayed in Figure 1. 

More general views of HA, including weapons and time periods, can be 
formulated as multicommodity compositions of generalized networks. 

The preceding perspective of the HA problem as a generalized network 
is of more than purely academic interest. Work has been published 
[Ref. 15] and research continues [Ref. 16] solving generalized networks 
with side constraints (such as flight composition constraints). A 
commercial quality optimization system exploiting this new technology 
could provide a new vehicle for enhancing the efficiency with which HA is 
processed . 
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AIRCRAFT/WEAPON 

i 



TARGET 

j 




ARC MULTIPLIER'. Pm 
ARC COST '.-V| ^ 



ARC COST: 0.0 
ARC MIN'. /) 
ARC MAX'. Tj 



Figure 1. Generalized Network View of the Simplified HA Internal Model 

B. AN EXPANSION OF THE HA INTERNAL MODEL TO INCLUDE WEAPONS AND TIME 
Headquarters Armament Division (AFSC/XR) proposed to Headquarters 
United States Air Force in January 1982 that HA be reformulated as a 
nonlinear mixed integer problem. Appendix G presents the basic concept 
for the reformulation. As reported by Dean [Ref. 17: pp. 56, 59], the 

X-System has solved nonlinear integer and mixed integer programs. 

(However, the size of the integer problems Dean reported were considerably 
smaller than the typical HA internal model.) The proposed formulation is 
basically an expansion of the dimension of the internal model to include 
weapons. 

Unfortunately, the model provided in Appendix G is intractable 
as stated and does not address all the issues pertinent to HA. It 
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entirely overlooks two model components described previously; weather 
and optimization over time. Also, it cannot ensure weapon selection of a 
"best" weapon for a particular aircraft/target combination prior to a 

a. U| 

"second best" weapon and so on to the "n best." A prioritization 
constraint is presented composed, in part, of binary indicator variables. 
However, the mechanism for enforcing the use of these binary variables is 
not presented. 

While a device for enforcing the use of the binary variables can be 
stated, this problem may be solved at less expense by allowing the 
selection of sorties to occur over the entire range of weapons considered 
by the problem. If the selection priority is ordered by (expected 

type j target kills by a type i aircraft loaded with type k weapon), 
prioritization is then enforced intrinsically by the presence of P... 

1 J K 

in the objective function. This alleviates the need for both prioritiza- 
tion constraints and binary variables, so that by expanding the size of 
the problem, the use of binary variables can be eliminated and the 
problem is made much easier to solve. 

Using the AFXR/SC proposal, without integer variables, as basis for a 
new approach, optimization over time could be addressed by further 
expanding the dimension of the problem to include time. The resulting 
formulation would be: 



Maximize: 

^ikn 



J 

E 

j=i 



E V. (K. ) 

jn' jn^ 



subject to; 
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1 . 



Sorties available. 



n NW / J n-1 

(1 - r) RS,^ < ^ E . E^ F(jkp Sijk. 



< (1 t r) E K5 

n=l 

2. Target constraints. 



In 



1, •••) I> n 1, 



N; 



-jn - ^j^^ijkn^ijkn - *^in’ 



jn’ 



1, ...J; n = 1, ..., N; 



3. Weapon logistics . 



I J N 

Z- S ^ik ^i jkn — '*^k ’ k = 1, ..., NW; 
i=l j=l n=l ^ 



4. Flight composition criteria 
5 



m 



, , y / , 



k = 1,2, ..., k; m = 1, 2, 



M; 



5. Target losses . 



K* “D*/ j_i\» 

jn j(n + 1) 

All parameters retain their previous definitions unless (re)defined 
in the following: 



T. -C. I NW 

^ " ill 1^1 ^ijkn^ijkn^^^’ 



jn 



k = weapon type index (k = 1, ..., NW); 

n = time period index (n = 1, ..., N); 

r = sortie equality range restriction factor; 

W|^ = number of available weapons of type k; 

= resupply quantity (initial quantity, if n = 1) of sorties 



available for aircraft i in period n; 
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jn 






L.|^ = number of weapons of type k loaded on aircraft of type i; 

*^ijkn " i^i^iendly aircraft attrition factor. 

Time periods would be connected by a friendly aircraft attrition 

factor, associated with each aircraft/target/weapon combination. 

The target losses constraints equate dead targets in period n + 1 (^j(n+i)) 

with those killed by the end of period n (K. ). 

jn' 

Unfortunately, if any C. ^ 0, the target losses constraint is 

J 

nonlinear and the model becomes somewhat more difficult to solve. 

Also, the attrition factor, might realistically be defined as a 

nonlinear function of time period or prior attrition. However, it is 
believed these nonlinearities can be accommodated by a state-of-the-art 
optimi zer . 

The Weapons Logistic constraint is appropriate when HA is used as a 
munitions procurement model. In HA's secondary role, providing theatre 
commanders with insights into scenarios, a weapon logistics resupply 
constraint might be preferred. Let Wj^^ represent the resupply of weapon 
type k arriving in the theatre at the beginning of time period n. Such a 
constraint can be stated as follows: 

3.1 Weapon Logistics resupply . 

I J 

n £ CN i=l j=l n £ CN ^ 

CN = 1, ..., N; k = 1, ..., NW. 

However, this constraint, if used as stated, produces a dense problem 
matrix. Model clarity and optimization performance would improve with 



36 



the addition of variables i 0 |^^ to represent the unused weapons of type k 
remaining at the end of time period n. This new constraint can be stated 
as follows: 

3.1.2 Weapon logistics resupply . 

I J 

“k(n-l) ’^kn " ‘’ik^ijkn " “kn’ 

n = 1, ..., N; k = 1, ..., NW 

where. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



This research was undertaken with the belief that HA was a single 
optimization model. Two models were identified and characterized. An 
internal (optimization) model was found, and an accompanying external 
model was discovered to principally constitute a circumvention of internal 
model shortcomings. Both models have been analyzed in detail, with special 
emphasis on the context of their use: scenario evaluation for theatres. 

The potential strengths and weaknesses of HA have been investigated. 

HA has been provided with a new, fast optimizer. A new prototype 
external model. Weather Wars, which corrects some of the more obvious 
shortcomings of HA, has been designed and coded. To a degree, the 
prototype suffers from the same lack of standard calling disciplines, 
variable naming conventions, and documentation as its predecessor, 
precisely because it is a modification of its predecessor. Any futher 
attempts to modify the existing external model are unlikely to provide 
continuing user satisfaction. However, an entirely new HA can be written 
to exploit the insights gained in this research. 

What HA has always needed is a single model, a single formulation, 
which encompasses enough of the real problem to be adjudged realistic. 
While this may have been impossible in 1974, given the state-of-the-art 
in optimization, it is feasible now. A single model such as that proposed 
in the preceding section would not only be easier to understand, but when 
implemented should actually run faster than the current HA, providing 
better solutions. 
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APPENDIX A. LINEAR RESTATEMENT OF TARGET CONSTRAINTS 



Target constraints are of the form: 



1 1 ’’’j* j = 1, 2, J; 



which is restated: 



j “ 1, 2, • • • 9 






First the upper inequality is solved for the linear term z 



-C 



I 



-C. I 

(1 - exp ( -^ (a +E P..S..))) < c., 
j ^ i = l “J 

1-C.<exp . E P,jS,j)». 



log (1 - Cj) < ( :^ . Ej 



■T. 



I 



-ji(log (1 - Cj)) > (aj . _E P,jS,j), 

TT (’09 <1 - ■ ’j - ,?1 

Similarly, the lower inequality is solved: 

T. -C. I 

S- U - exp ( ^ (a + 21 PiiSii))), 



i. yi < (1 - exp ( ^ (a + Y. PiiSii))), 
J 'j 'j J i=i 
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C. -C. I 

1 - ^ > (exp ( -^ (a + £ p. .s..))), 

J ij - ij J iJ iJ 

C. -C. I 

log (1 . i J) > ( -^ (a + £ J), 

J Ij Ij J 1 = 1 U 1J 



-T. 

"C^ (log (1 
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APPENDIX 8. ALGEBRAIC MANIPULATION OF K. - 0- 

J J 



If 



T. -C. I 

C = J (1 - exp ( -yi (a . + £ Pi -S-.))), 

J Cj Tj J U 



then noting 



-C . C . 

exp ( ^ a ) = 1 - yi D . , 

j j 



and letting 



X, = E P.,S,,, 



i=l 






the term K. - D. (from the current HA internal model objective function) 
0 J 

can be restated as follows: 

- “j 'tf ‘ - “j 

= i(l-(l-JDj)(exp(:^iXj))).Dj(^)(i) 



T 

_ J 



C -C c 

J (1 - (1 - j 0.)(exp ( ^ X,)) - 0,( J- )) 
Pj ij J ij 0 0 ij 



T C C -C 

= j (d - 0 .fl ) - (1 - yi D )(exp (X, )) 

j J 0 ^ ^ J 

T. C. T. -C. 

° ■ “j >> 



= (X0j)-(^-O.)(exp(Xj:^)) 

J J J 
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T. 

= ( 7^ - D.)(l - exp 





)). 



If C. has a value of one, then K. - D- simplifies further to: 

K. -D^ = (Tj-Dj)(l-exp-i-)). 

If Cj has a value approaching zero, then Kj - approaches Xj. 
Recalling (1 - exp (-x)) is approximately x for small x, and noting 
-C. 

( X^) -> 0 as Cj -> 0, 

and further noting that as C. approaches zero, the quotient T,-/C- 
will become so large as to leave the effect of D. negligible, K- - 0- 

\J \j \j 

then becomes 



J J 

In a target rich environment, no matter what the value of C., X./T. 

\J \J \J 

will be quite small. Noting the fact that (1 - exp (-x)) is approximately 
X for small x, K. - D- is finally approximated by: 

X. X.D. 

- “:•> tJ = 



X. 



= Xj, for small yi . 



Thus, Kj is nearly linear in a target rich environment. 
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APPENDIX C 



ORIGINAL HEAVY ATTACK INPUT-PROCESS-OUTPUT 




MUNITIONS MIX 



TARGETS KILLED 
IN interval r 



K£ Y: ^ A Input/ Oufpwf 

^ -y factor 



Processing 

function 
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APPENDIX 0. 



HIERARCHY CHART OF HA PROGRAMS IN USE AS OF 1 JUNE 1982 
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APPENDIX E. VARIABLE NAME CHANGES IN HA PRIOR TO OPTIMIZATION 



Formul at ion 
Variable 


Subroutine 

NONLIN 


Subroutine 

OPTMUM 


Main 

Program 


Si 


SMAX 


S 


SORT 


I 


XNAT 


NAC 


NAC 


J 


XNTT 


NTGTS 


NTGTS 


b 


T 


TTGT 


— 


“j 


D 


DTGT 


TD 


l . 
J 


SL 


TL 


— 




CC 


CC 


CC 


''j 


VT 


VT 


VT 




PROS 


EKS 


EKS 


no. side 
constraints 


XNADC 


NADC 


— 


m 


ROW 


IDAC 


— 




PCTX 


PCT 


— 


# in set J 

m 


NTSX 


NUMTGT 


— 




lARRY 


lARRY 


— 


'ij 


YS 


SOL 


MSORT/SOL 


Obj Fen Value 


GM 


OBJ 


OBJ 
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APPENDIX F. NAP MODELS INPUT-PROCESS-OUTPUT RELATIONSHIPS 



File 1 

Weapon delivery profile data 



File 2 

Weapon effectiveness 



File 3 
Payloads 



File 7 

Fractional kills per pass (FKP) 



File 5 

Target description 



! 

File 4 

Fractional kills per pass (FKP) 
by weapon/delivery condition 
for each aircraft and target 
in-year weapon 





File 3 






Payloads 





Y 

File 7 

FKPs for each aircraft and 
weapon for each target 



File 8 

Attrition for each profile and 
aircraft against every target 



File 10 

Weapon profile and bomb damage 
assessment data 



File 4 

Acquisition probabilities 






WEAPONICAL 

Computes the fractional kills 
per pass (FKP) for every 
aircraft delivery condition 
and target combination. 



I INBETWEEN 

I Reformates FKP data into one 
file with the FKPs grouped by 
1 targets 



SURVIVOR 

Computes the expected kills 
per sortie (EKS) and the 
attrition sustained by attack- 
ing each target by a specified 
weapon and delivery profile 
with a specified aircraft 
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T 






le 2 

EKS and attrition for each 
aircraft and weapon by target 



File 5 

Cost data. Weather states. 
Delivery conditions. 



I 

I 

File 8 

Preferred weapons by delivery 
band for each target 



File 5 

Number of sorties. Target 
values. Constraints. 







SELECTOR 

Selects the preferred delivery 
condition/weapon for each air- 
craft/target/del ivery band 
combination. 



i 

T 

I Rank ordering 
I of weapons 



I HEAVY ATTACK 

I Identifies the sortie alloca- 
I tion (aircraft, weapons, and 
I targets) that maximizes 
I mi 1 i tary worth 



I Sortie, target and aircraft 
1 allocations that maximize 
I military worth 
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APPENDIX G. AFSC/XR PROPOSED NEW HA METHODOLOGY 



Definition of Variables: 



i 



j 

k 






^ijk 



V. 

J 






K 



L 



ijk 



n 



n 



(i J) 



P 



m 



i 



m 



Aircraft type index (i =1, I); 

Target type index (j = 1, J); 

Weapon type index (k = 1, K); 

Number of j type targets; 

Cumulative number of targets killed in prior time periods; 

Expected number of type j targets killed per type i aircraft 
loaded with type k weapon when no other targets of type j have 
been previously killed and when conditions of confirmability 
are perfect; 

Number of sorties of type i aircraft flown with weapon 
type k against target type j; 

Quantity of sorties (fragable) available for aircraft type i; 
Value (Military Worth) of target type j; 

Lower bound on targets of type j to be killed; 

Number of type j targets killed; 

Target kill conf irmabi 1 ity parameter for j type targets 
(as defined in old formulation); 

Quantity of type k weapons available; 

Standard loadout of weapon type k on aircraft type i used 
against target type j; 

Number of members in the weapon prioritization sets; 

The p priority member (a k value) of the ordered weapon 
priority set of k values valid for the i^” aircraft type, 
j target type combination (p ranges from 1 to m^^ best); 

flight composition constraint index (m = 1, 2, M); 

aircraft type; 
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5^ = +1, maximum {-1, minimum); 

= set of targets for which a maximum (minimum) flight 
composition is required; and 

9^ = maximum (minimum) proportion of sorties flown by aircraft 

type i^ against targets included in set 

The formulation: 



J 

Maximize Z = 

j=l ^ ^ ^ 



where 



a . = 






/ 


-C. 


(l - exp 


J 

T. 


\ 


J 




C.D. 


log (1 - 


J J 
“T7 






I 

E 

i = l 






subject to: 

Sortie constraint; 

J K 

L E < F.; i = 1, I; 

j=l k=l ^ 

Target constraints; 

1 Q j 1 j 5 j = 1 , . . . , J ; 

Side (flight composition) constraints; 



m 



Ed 

jeJ 



m 






ik — 

m = 1, 2, 
k = 1, 2, 



Weapon constraints; 



I J 

^E ^E 1 k = l, ...,K; 



M; 

K; 
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Prioritization constraints; 






1 = 1 , 
j = 1, 
P = 1, 



where this indicator variable is defined as: 






P+1 



0’ ^iin 
1; otherwise 



= 0 



I; 

J; 

n; 
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